is possible to account for all of these constraints in a self-consistent manner. 
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The wealth of observational data about Jupiter and Saturn provides strong 
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c/3 . 1. Introduction 

.^J [ The standard giant planet formation scenario is the so-called core-accretion model. In 

this model, a solid core is formed first by the accretion of solid planetesimals. As the core 
grows, it eventually becomes massive enough to gravitationally bind some nebular gas, form- 
ing a gaseous envelope in hydrostatic equilibrium. The further increase in core and envelope 
masses lead to larger and larger radiative losses which ultimately prevent the existence of 
an equilibrium envelope. Runaway gas accretion occurs, rapidly building up a giant planet. 
This scenario, which naturally implies giant planets enriched in heavy elements compared to 
the Sun, has suffered so far from the problem that the resulting formation time is comparable 
to, or longer than, the lifetime of protoplanetary disks as inferred from observations (Pollack 
et al. 1996 - hereafter P96, Haisch et al. 2001). One approach to solve this long-standing 
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problem has been to revise the opacities used to model the planet's envelope (see Hubickyj 
et al. 2003) and/or to allow for a local enhancement in the number of planetesimals (Klahr 
& Bodenheimer 2003). However, Alibert et al. (2004,2005a - hereafter A04 and A05a) have 
shown recently that extending the original core accretion scenario to include migration of 
the growing planet and protoplanetary disk evolution, results in a formation speed-up by 
over an order of magnitude even without local density enhancements or modified opacities 1 . 
In this Letter, we show that in addition to solving the formation timescale problem, these 
models can also account for the characteristics of the two most well known giant planets: 
Jupiter and Saturn. 

Within the framework of our model, we find that the uncertainties on the characteris- 
tics of the initial protoplanetary disk are large enough to allow us to match the observed 
properties of a single planet (Jupiter for example) relatively easily. The situation is more 
complicated with two planets forming within the same protoplanetary disk. For each satis- 
factory model matching the observed properties of Jupiter (total mass, distance to the sun, 
mass of the core, total mass of heavy elements and volatiles enrichments), only two parame- 
ters are left in order to account for the same five quantities in Saturn: the initial location of 
the embryo and the time offset (that can be equal to 0) between the start of the formation of 
both planets. The purpose of this Letter is to show that, by assuming reasonable properties 
for the initial protoplanetary disk, it is possible to construct models of Jupiter and Saturn, 
compatible with all the observations detailled in the next two paragraphs. 

Using measurements of Jupiter and Saturn (mass, radius, surface temperature, gravi- 
tational moments, etc.) and state-of-the-art structure modeling, Saumon & Guillot (2004 - 
hereafter SG04) have derived important constraints regarding the possible internal structure 
of these planets. From this modeling, M core , the mass of the core of the planet and M Zienve , 
the amount of heavy elements in the envelope (assumed to be homogeneously distributed) 
can be obtained. In the case of Jupiter, the maximum total amount of heavy elements 
present in the planet (M corc + Mz jenvc ) is of the order of ~ 42M e (Earth masses), whereas 
the mass of the core can vary from to 13M e . This large uncertainty is essentially due to 
the undertainties in the equation of state of hydrogen. In the case of Saturn, Mz iCnvc ranges 
from nearly to 10 M®, the mass of the core being between 8 and 25 M®. Note however that 
the mass of the solid core might be reduced by up to ~ 7M e depending upon the extend of 
sedimented helium, a process which is required to explain the present day luminosity of the 
planet (Fortney & Hubbard 2003, Guillot 2005). 



1 The capability of migration to prevent the isolation of the protoplanetary core was already stated by 
Ward (1989) and Ward & Hahn (1995). 
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Abundances of some volatile species in the atmosphere of Jupiter have been measured 
using the mass spectrometer on-board the Galileo probe (see Mahaffy et al. 2000, Wong 
et al. 2004). These measurements show that the planet's atmosphere is enriched in C, N, 
S, Ar, Kr, and Xe by a factor of 3.7 ± 0.9, 3.2 ± 1.2, 2.7 ± 0.6, 1.8 ± 0.4, 2.4 ± 0.4, and 
2.1 ± 0.4 respectively compared to solar values (Lodders 2003). For Saturn, ground-based 
observations (Brigg & Sackett 1989, Kerola et al. 1997) have shown that C and N are 
enriched by a factor of 3.2 ± 0.8 and 2.4 ± 0.5 compared to the solar values. Since the two 
planets are almost entirely convective, we assume that these enrichments are representative 
of the mean envelope composition. 

In Sect. 2 of this Letter, we give a short presentation of our formation models. In Sect. 
3, we apply these models to Jupiter and Saturn, and show an important role of Jupiter's 
formation on that of Saturn. In Sect. 4, we calculate the enrichments in volatile species in 
the atmosphere of the two planets, and Sect. 5 is devoted to summary and conclusions. 

2. Formation models 

Our formation models consist in the simulation of the time evolution of the proto- 
planetary disk and of the two planetary seed embryos that will eventually lead to Jupiter 
and Saturn. We calculate in a consistent way the structure and evolution of the disk, the 
migration of the planets, and their growth in mass due to accretion of gas and planetesimals. 

The evolution of the protoplanetary disk is calculated in the framework of the a formal- 
ism (Shakura & Sunyaev 1973). The initial gas surface density £ inside the protoplanetary 
disk (which extends from 0.25 AU to 30 AU) is given by S oc r~ 3 / 2 . The gas to solids 
ratio is constant in the whole disk (the embryos always stay beyond the ice line), with 
E gas /E solids = 70. The gas surface density evolves as a result of viscous transport and pho- 
toevaporation: 

+ t w {r). 

The photoevaporation term S w (r) is taken as in Veras & Armitage (2004). The thermo- 
dynamical properties of the disk as function of position and surface density (temperature, 
pressure, density scale height), as well as the mean viscosity u, are calculated by solving the 
vertical structure equations using the method presented in Papaloizou & Terquem (1999) 
and A05a. These quantities are used to determine the composition of the ices incorporated 
in the planetesimals, and finally the enrichments in volatile species in the two planets (see 
Sect. 4). 
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The key point in our models is that both planets are formed concurrently within the 
same disk, hence the physical assumptions and initial properties of the nebula are the same 
for both. For a given disk model, we begin by searching a satisfactory model matching 
the observed properties of Jupiter. Once such a model is found, we try to adjust the two 
remaining parameters (initial location of the embryo and time delay) to find a similarly 
suitable model for Saturn. Our entire approach, as well as some tests we have made to 
check our code, can be found elsewhere (A05a), we give here some details on two points, the 
calculation of M corc and M Z en ve, and the migration rates. 

To estimate M core and Mz ien ve, we compute the fate of the infalling planetesimals by 
computing their trajectory inside the envelope as well as their mass loss. The latter results 
from thermal effects as well as mechanical ablation due to Rayleigh- Taylor instabilities on 
the planetesimals' surface (Korycansky et al. 2002). This allows us to determine the fraction 
of planetesimals' mass that directly reaches the core, which we identify as M core , the core 
mass at the end of the formation process. The mass deposited inside the envelope is assumed, 
due to convection 2 , to be homogeneously distributed within the envelope and is identified as 
Mz^enve, the mass of heavy elements in the planet's envelope. Finally, we note that processes 
like core erosion or settling could occur during subsequent evolution of the planet and modify 
significantly the values of M core and M Z enve found here (see SG04). On the other hand, the 
total mass of heavy elements in the planets should remain constant. 

The calculation of the planetesimal's trajectory also gives the place where the energy of 
planetesimals is deposited. This quantity is used to calculate the structure of the forming 
planet, by solving the standard internal structure equations: the amount of energy released 
by infalling planetesimals enters in the energy equation, in the sinking approximation (see 
P96 and A05a). 

Low mass planets undergo type I migration at a rate being linear with the planet's 
mass. However, the most recent analytical estimates of type I migration rates by Tanaka 
et al. (2002), which have been derived assuming a laminar disk, are much too large to be 
compatible with the observed frequency of extra-solar planets. Therefore, planet survival 
implies a significantly reduced rate of type I migration. First hints how this could be achieved 
have been obtained by Nelson & Papaloizou (2004) in numerical modelling of turbulent disks 
in which much reduced migration rates have been obtained. In our calculations, we have 
reduced the rate of type I migration by multiplying the analytical estimates by an arbitrary 
factor fi, whose value varied in order to check its influence on the results. For higher mass 



2 The interior of the planets is largely unstable with regard to convection, even with the presence of 
molecular weight gradients: the radiative gradient is dominant over the adiabatic one by orders of magnitude. 
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planets, the migration is of type II (Ward 1997), the rate being independant of the planet 
mass. When the mass of the planet becomes comparable to the one of the disk, migration 
slows down and eventually stoppes. The switch from type I to type II occurs when the Hill's 
radius of the planet is equal to the disk density scale height, which is calculated with the 
vertical structure of the disk. Finally, note that we do not take into account gravitational 
interactions between the two forming planets that could alter the migration rates. 

3. Jupiter and Saturn formation 

We consider values of fi between (no type I migration) and 0.03 (as we shall see below, 
higher values would imply too large starting locations of proto- Jupiter to account for the 
present structure of Saturn). For this range, we find suitable Jupiters to form from embryos 
starting between 9.2 AU (Astronomical Units) and 13.5 AU in a disk with a total mass 
ranging from 0.05 to 0.035 M (solar masses) and a total photoevaporation rate comprised 
between 1 and 1.5 xlO~ 8 M / yr. For all the cases considered here (/i=0, 0.001, 0.005, 0.01 
and 0.03), it was possible to form within 3 Myr a planet whose final mass, location and global 
internal structure were compatible with the Jupiter ones (see Fig. Id). We note that the 
final structure of the planet (M core , Mz envc ) is independant of the assumed type I migration 
rate. This rate only gives the starting location of the embryo. 

Concurrently with the formation of Jupiter, we also follow the growth of the proto- 
Saturn embryo. The latter is started at a larger heliocentric distance and with an arbitrary 
time delay. This implies that, depending upon initial conditions, proto-Saturn may actually 
enter a region of the disk already visited and consequently modified (less planetesimals for 
example) by proto- Jupiter (see Fig. lc). In Fig. 2, we present the successful Saturn forma- 
tion model corresponding to the Jupiter model presented in Fig. 1 (red curves, /i=0.001). 
The synthetic Saturn started as an embryo at 11.9 AU, 0.2 Myr after proto- Jupiter. The 
resulting planet exhibits characteristics quite comparable to the actual Saturn (see Fig. 2d). 
The mass of the core is slightly lower than the one allowed by SG04. However, we recall 
that the mass derived in SG04 may be decreased by up to ~ 7M e due to the sedimentation 
of helium (see Sect. 1). 

The mass of Saturn's final core is similar to the one obtained for Jupiter. This is because 
the core is built from the infalling planetesimals that are able to traverse the gaseous envelope 
without being disrupted. For a fixed envelope, disruption is essentially a function of the size 
of the planetesimals which in our work is assumed to be identical at all locations (100 km). 
Increasing the mass of the planetesimals by a factor ten leads to core masses of the order of 
8M ffi . 
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The importance of Jupiter's wake on the formation of Saturn is evidenced by the green 
curves in Fig. 2a which were obtained by forming Saturn in the absence of Jupiter in an 
otherwise identical disk. In this case, the increased rate of planetesimal infall prevented the 
accretion of a sizeable envelope, and the resulting planet, while at Saturn's current location, 
remained quite small (20M©). This can be explained by the fact that the gas accretion rate is 
inversely dependent upon the energy deposited by infalling planetesimals. Thus, once proto- 
Saturn's feeding zone enters a region previously depleted in planetesimals by the passage of 
Jupiter, their infall is reduced and gas accretion proceeds at a faster rate, ultimately leading 
to a more massive planet than in the case without planetesimal depletion. In the latter 
case (ignoring the effects of Jupiter's formation), we checked that even by varying the initial 
location and formation starting time, it was not possible to obtain a Saturn-mass planet at 
its current location (see Fig. 2c). Increasing the starting location of Jupiter (beyond ~ 10 
AU, corresponding to f\ larger than 0.01), results in Saturn-like planets containing too few 
heavy elements compared to the actual planet. Moreover, the mass of accreted planetesimals 
never reaches a level which could trigger a significant accretion of gas. At the present location 
of Saturn, the synthetic planet remains less massive than the actual one (see Fig 2c, blue 
curve). 



4. Enrichments in volatile species 

We now concentrate on the models that can form both Jupiter and Saturn (red ones 
in Fig. la and Fig. 2a) and examine whether they can also account for the volatile abun- 
dances measured in the atmospheres of the two planets. To do this, we use the clathrate 
trapping theory (Lunine & Stevenson 1985) and the thermodynamical conditions inside the 
disk as calculated in our models to compute the composition of ices incorporated in the 
planetesimals. Knowing the mass of accreted planetesimals, we compute the total expected 
abundances of some volatile species, in our case C, N, S, Ar, Kr, and Xe. 

We use for this calculation the recent solar abundances determinations of Lodders 
(2003). C is set to have been present in the solar nebula vapor phase under the form of 
CO2, CO and CH 4 , with CO2:CO:CH 4 =30:10:l, ratios which are compatible with ISM mea- 
surements (see Allamandola et al. 1999, Gibb et al. 2004). Moreover, N is taken to have 
been present under the form of N 2 and NH 3 , with a ratio NH 3 :N 2 =1, and S under the form 
of H 2 S and other sulfur compounds (Pasek et al. 2005). Other initial ratios of C02:CO:CH 4 
and NH 3 :N 2 can lead to slightly different abundances of volatiles, but do not modify our 
main conclusions. We note finally that C0 2 crystallizes as a pure condensate prior to be 
clathrated (see Alibert et al. 2005b) which has a considerable influence on the total amount 
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of water required for trapping all the volatiles in the planet. 

The results of these abundances calculations in the case of Jupiter have been presented 
elsewhere in details (see Alibert et al. 2005b), we only summarize the main conclusions 
here: C, N, S, Ar, Kr and Xe are enriched respectively by a factor of about 2.8, 2.5, 2.1, 
2, 2.1, 2.6 compared to their solar values. These values are compatible with the in situ 
measurements made by the Galileo probe and recalled at the beginning of this Letter. The 
resulting enrichment for oxygen (not yet measured) is at least O/H = 3.4 x 10~ 3 or ~ 6 
times the solar value. 

In the case of the Saturn formation model, we obtain enrichments of 2.4 and 2.2 com- 
pared to solar values for C and N respectively. This is again compatible with the ground- 
based observations quoted in the introduction. Moreover, we predict that S, Ar, Kr, and 
Xe are enhanced by a factor of respectively 1.9, 1.7, 1.9 and 2.3. Our formation model pre- 
dicts the accretion of ~ 13.2M e of heavy elements, and the trapping of the volatiles result 
from the accretion of at least 5.4M® of ices (depending on the efficacity of the clathration 
process). These two calculations imply that the mean Ices/Rocks (I/R) ratio of accreted 
planetesimals was > 0.7, a value consistent with the one inferred for Saturnian satellites. 
Finally, the resulting enrichment of O in Saturn is O/H ~ 3 x 10~ 3 , ie 5.2 times the solar 
value. 



5. Summary and discussion 

We have calculated in this Letter formation models of the two gas giant planets of our 
Solar System, in the framework of our extended core-accretion models taking into account 
migration and disk evolution. The calculations presented here are simplified in some aspects, 
that could be improved in the future. In particular, our disk model is calculated in the 
framework of the a formalism of Shakura & Sunyaev (1973), which itself is a limitation. 
Moreover, we do not take into account gravitational interactions between the two planets 
that can alter the migration rates. 

Our calculations allowed us to give an estimate of the core mass, enrichment in heavy 
elements, and enrichment in volatile species that can be compared with observational data 
about Jupiter and Saturn. These calculations therefore show that our models can lead to the 
formation of two giant planets closely resembling our Jupiter and Saturn in less than 3 Myr. 
In order for our synthetic planets to match the bulk properties of the two gas giants in our 
solar system, we found that Jupiter must have started at a heliocentric distance smaller than 
~ 10 AU otherwise Saturn, which follows in its trail, cannot accrete enough heavy elements. 



- 8- 



The heavy elements content as well as the core mass obtained for our synthetic planets are 
in good agreement with the interior models of SG04. 

Finally, in the framework of our model, the enrichments (compared to solar) of some 
volatile species measured in the atmosphere of both giant planets can also be accounted 
for in a self-consistent manner. However, we note that, recently, the Cassini spacecraft 
measurements have led to a revised value of the abundance of C in Saturn's atmosphere 
of 8.1 ± 1.6 times the solar value (see Flasar et al. 2005). Using the clathrate hydrate 
trapping theory, and assuming the minimum value for the C abundance (6.5 times the solar 
value), we obtain abundances of N, S, Ar, Kr, and Xe of respectively 5.9, 5, 4.7, 5.1 and 
6.1 times the solar values. The resulting O enrichment would be about 14 times the solar 
value. These predicted enrichments are significantly higher than the ones quoted in Sect. 
1, in particular for N. The future confirmation of both the new measurement of C and of 
the "old" value of N would imply that there has been some unknown fractionation processes 
between these species in the solar nebula gas-phase, and that their abundances in Saturn 
cannot be explained using solely the standard clathrate hydrate trapping theory. However, 
the measurement of C in Saturn by the Cassini spacecraft may be subject to revisions in a 
near future. 

This work was supported in part by the Swiss National Science Foundation. O.M. was 
supported by an ESA external fellowship. 
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Fig. 1. — Jupiter formation models. The red (resp. blue) curves are obtained using j\ = 
0.001 (resp. 0.03). (a) Mass of accreted planetesimals (dashed lines), total mass (dotted 
lines) and mass of solid core (solid lines) as a function of time for two simulations, (b) 
Heliocentric distance as a function of time, (c) Final surface density of planetesimals S p 
as a function of heliocentric distance, (d) Core mass (M core ) and mass of heavy elements 
(Mz,enve) dissolved in the envelope. The black curve gives the domain allowed by the present 
day structure models of SG04. If some core dissolution occurs after the formation, the two 
points would evolve along the lines. For all the values of fi we used, the points representing 
the final structures are very close to the red and blue ones. These points are not represented 
here for clarity. 
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Fig. 2. — Saturn formation models. The red (resp. green) lines are obtained while taking 
into account the effect of Jupiter (resp. without this effect). (a,b) Same as Fig. 1 a,b. (c) 
Red small dots: different Saturn models (final mass and position), varying the initial position 
and the time delay after the beginning of Jupiter's formation. The approximate maximum 
mass that can be reached by proto-Saturn at a given distance to the sun is indicated by 
the lines. The blue line is similar to the red one, except that the corresponding Jupiter 
formation is calculated with fj = 0.03 (blue curves in Fig. 1). The green line is similar to 
the red one, except that the effect of Jupiter's wake on Saturn formation is not taken into 
account. The black star indicates Saturn's position in this diagram, (d) Core mass (M core ) 
and mass of heavy elements (M Z enve ) dissolved in the envelope. The black curve gives the 
domain allowed by the present day structure models of SG04, taking into account a possible 
shell of sedimented helium, of mass between and 7M ffi . 



